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We present the magnetic phase diagram of artificial H2 molecule in lateral quantum dots doped 
with magnetic impurities as a function of external magnetic field and plunger gate voltage. The 
onset of Mn-Mn antiferromagnetic-ferromagnetic transition follows the electron spin singlet-triplet 
transition. We deploy a configuration-interaction method to exactly diagonalize the electron-Mn 
Hamiltonian and map it to an effective Mn-Mn Heisenberg Hamiltonian. We find that Mn-Mn ex- 
change coupling can be described by RKKY-interaction/magnetic-polaron in weak/strong electron- 
Mn coupling at low/high magnetic fields. 

PACS numbers: 75.75.+a,75.50.Pp,85.75.-d 



I. INTRODUCTION 



There is cu rrently significa nt experimenta l 1 * 2 * 3 * 4 * 5 * 6 * 7 ! 
and theoretica l 8 * 9 * 10 * 11 ** 13 * 14 ! interest in semiconductor 
quantum dots (QDs) doped with magnetic impurities. 
Fabrication of hybrid systems consisting magnetic ions in 
a controlled electronic environment provides an interest- 
ing interplay between interaction effects and magnetism. 
In particular the application of spin of electrons in quan- 
tum dot molecules for generation of electron entangle- 
ment and quantum information processing in solid state 
devices is of current interestP^EHl One of the challenges 
in use of spin of electrons in scalable quantum computer 
devices is the spin dephasing due to interference by spin- 
orbit coupling and nuclear hyperfine interaction^. The 
broadening of the electron envelop wavefunction in a QD 
is determined by the length scale of the confining poten- 
tial that is comparable with the size of QD (typically 
« 1-1000 nm). In materials with high abundance of 
spin-carrying nuclear isotopes, the electron interacts with 
large number of nuclei in host semiconductor. Although 
the strength of nuclear hyperfine coupling is small (typ- 
ically ss 1/ieV) compare to other relevant energy scales, 
but because of the broadening of the confined electron in 
nano-meter length scale, it interacts with a large number 
of nuclei that contribute to the electron spin dephasing. 
The use of magnetic moment of nuclear impurities in 
host semiconductors for the quantum information pro- 
cessing have appeared to overcome this limitation 1 18 * 19 !. 
Recently the system of 13 C atoms in two-electron nan- 
otube quantum dot molecules has been studied 20 . The 
advantage of using singlet eigen-states of coupled mag- 
netic moments of nucleus of 13 C in organic molecules 
and their effects in dramatic enhancement of the spin 
life-times needed for imaging the metabolic pathways in 
living systems by hyperpolarization methods has been 
investigated^. Moreover, molecules doped with 13 C has 
been used in demonstration of quantum teleportation us- 
ing NMR techniques^. Similarly the magnetic dipolc 
moment of nucleus of 3 ^-impurities in Si-based quan- 
tum computer model proposed by Kane^l appears to be 



promising in making quantum computer solid state de- 
vices. Similar to magnetic dipole moment of nucleus of 
13 C and/or 31 P, the magnetic moment of electrons in 
d-shell of magnetic impurities (such as Mn, Fe, Co) in 
lateral quantum dot molecules can be used for quantum 
information processing as they have been used for fabri- 
cation of molecular magnets^. 

In this work we focus on theoretical study of the phase 
diagram and spin transitions of coupled magnetic im- 
purities (e.g. Mn) doped in two electron quantum dot 
molecule, an artificial H 2 molecule. Here we calculate 
the Mn-Mn effective Hamiltonian mediated by electrons 
and show that the onset of spin-polarized state with 
finite Mn-magnetization, corresponds to spin singlet- 
triplet transition of two electrons in QD molecule. This 
transition is analogous to para-ortho transition of nu- 
cleus of solid H2 where the electron-nuclei hyperfine in- 
teraction opens the energy gap between para and ortho 
states of H 2 molecule 1 24 * 25 ! In the small (large) magnetic 
fields the spin singlet (triplet) state is the ground state of 
electrons^ and thus the ground state of the coupled Mn 
is described by hydrogen molecule para (ortho) state. In 
the small magnetic fields the Mn-Mn interaction induced 
by electrons is calculated perturbatively in the electron- 
Mn weak interacting limit. It can be described effectively 
by RKKY-coupling 23 . In large magnetic field the ground 
state of electrons is spin-triplet. The electron-Mn inter- 
action is strong and magnetic-polaron state form. In this 
limit the ortho-state of artificial H 2 molecule is stable. 
We show the Mn-Mn exchange coupling can be controlled 
by inter and intra dot correlations, external magnetic 
field and gate voltage. The dependence of Mn-Mn inter- 
action to the external gate voltage and magnetic field me- 
diated by spin singlet-triplet transition among electrons 
in QDs opens up the possibility in using Mn-magnetic 
moment as qubit in quantum computation purposes. In 
contrast to electrons confined in QDs, because of highly 
localized d-clcctrons, Mn's interact directly with signif- 
icantly smaller number of nuclei in host semiconductor. 
The Mn d-electrons also do not interact directly with 
the host semiconductor Rashba and Drcssclhaus spin- 
orbit coupling^, and thus their spin coherence life-time 
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is expected to be longer than the QD electrons. How- 
ever, Mn-Mn interaction mediated by electrons in QDs 
are still vulnerable to the electron spin dcphasing medi- 
ated by QD electrons due to their hyperfine and spin- 
orbit couplings. In analogous to the 31 P in Si system^, 
the spin coherence time in the Mn system is expected to 
be longer than the QD electron system. Further inves- 
tigation is required to make a quantitative dependence 
of spin decoherence time of Mn on the electron-nuclear 
hyperfine interaction and their spin-orbit coupling. 



II. HAMILTONIAN 

We represent magnetic QD molecule by the Hamil- 
tonian H = H e + H em + H m , describing contributions 
of interacting electrons, electron-Mn (e-Mn) exchange, 
and direct Mn-Mn antiferromagnetic (AFM) coupling, 
respectively. Electrons confined in quasi-two-dimensional 
quantum dots in a uniform perpendicular magnetic field 
can be described by the effective mass Hamiltonian 
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where 



T=^(^ + - c A(r ) ) + V(x,y) (2) 

is the single electron Hamiltonian in magnetic field. Here 
(f) = (x, y) describes electron position, V{r) is the quan- 
tum dots confining potential, A(f) = ^B x r is the vec- 
tor potential, and B is the external magnetic field per- 
pendicular to the plane of confining potential, m* is 
the conduction-electron effective mass, — e is the elec- 
tron charge, and e is the host semiconductor dielectric 
constant. Eiz = \g e ^B^izB is the Zeeman spin split- 
ting, g e is the electron g-factor in host semiconductor, 
\xb is the Bohr magneton, and a is the Pauli matrix. The 
single particle eigenvalues (e aa ) and eigenvectors (ip aa ) 
are calculated by discretizing T in real space, and diag- 
onalizing the r esulting matrix using conjugated gradient 
algorithir i 16 l 26 l The single-particle (SP) states can be 
used as a basis in configuration-interaction (CI) calcu- 
lation that allows to diagonalize Hamiltonian H. The 
details of CI method can be found in RefP3. Denoting 
the creation (annihilation) operators for electron in non- 
interacting SP state \a<j) by cj^ (c QCT ), the Hamiltonian 
of an interacting system in second quantization can be 
written as 
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where the first term represents the single par- 
ticle Hamiltonian Eq.(j2|, and V a a;0a-' ,/j.cr' ,va = 



I dfj dr'cp^^cp^ir^^^ip^ir'^air), is 

the two-body Coulomb matrix clement. We de- 
scribe e-Mn exchange interaction by H em = 
-JsdJ2ijSi ' ~ R-/)> where J sd is the exchange 

coupling between electron spin s*j at and impurity 
spin Mi at R j 10 * 11 ! Note that H em is analogous of the 
isotropic (contact) part of electron-nucleus hyperfine 
interactiorP-2, responsible for para-ortho energy gap of 
solid H2 (for comparison see for example eq. 121.9, page 
498 of Ref P^j) . In second quantization it can be written 
as 

Hem = E E ^^[JM«£r«*t " 
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where J Q/3 (R/) = J s df* a (R/)^(R/). Finally we de- 
scribe Mn-Mn direct exchange interaction and Mn- 
Zeeman coupling by H m = J2i r JfJ'Mi ■ Mjr + 
E/ g-mfJ-BMizB, where jfj, is the direct Mn-Mn AFM 
coupling^, resembling the direct dipole-dipole interac- 
tion, and g m is the Mn g-factor. 



A. Confining potential 

For numerical calculation we model quantum 
dots molecules by the following confining potential 

V(x,y) = V L exp[- (*+g +ya ] + V R exp[- ^ W ] + 

V p exp[— S3 - ]• Here Vl,Vr describe the depth 

of the left and right quantum dot minima located at 
x = —a, y = and x = +0, y = 0, and V p is the plunger 
gate potential controlled by the central gate^. For 
identical dots, Vl = Vr = Vq, and confining potential 
exhibits inversion symmetry. Our numerical results are 
calculated for parameters based on (Cd,Mn)Te QDs with 



J., 



0.015 eV nnr\ m* 



0.106, 



10.6, 



g e = —1.67 and the effective Bohr radius 
and Rydberg energy Ry* — 12.1 
the confining potential by Vq = - 



g m = 2.02, 
= 5.29 nm 
meV. We parametrize 
■10, a = 2, A = 2.5, and 



A p x = 0.3, Ap y — 2.5, in effective atomic units. V t 



p ■ 



which controls the potential barrier, varies to control 
the inter-dot correlations, independent of the locations 
of the quantum dots. The choice of parameters ensures 
weakly coupled quantum dots. 



III. TWO LEVEL SYSTEM 

For the purpose of this work we consider a coupled 
quantum dot system filled with two electrons. It is con- 
venient to project the Hilbert space of two electrons into 
a two level system. The construction of two level system 
based on single particle orbitals localized in each dot is 
made by defining a pair of bonding-anti-bonding single 
particle orbitals <p±(r) = [ipt{T) ± ^fl(r)] /\/2(l ± W), 
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FIG. 1: (Color online) (a) Shown singlet-triplet energy gap 
(Ajy) of two electrons in a coupled quantum dot molecule as 
a function of plunger gate voltage (V p ) at B — 0. (b) Shown 
a 2 as a function of V v . 



FIG. 2: (Color online) Spin singlet-triplet phase diagram cal- 
culated by configuration interaction method for two electrons 
in lateral quantum dot molecules. 



where <pL{R){?) is the spatial part of SP wave-function 
localized in L (R) dot, and W — Re((L\R)) is the 
overlap integral. At zero (finite) magnetic field the SP 
orbitals are real (complex) functions. At H em = B = 
the lowest energy many body wave function (ground 
state) of two electrons is spatially symmetric with parity 
+ 1. Thus spin state of the ground state must be singlet: 
*G(ri,r 2 ) = [a(p + (f 1 )ip + (f 2 ) + /3<£_(ri)^_(f 2 )] \S ), 



where = \f\ 



and |5 ) = (| U) - I IT»A/2 



corresponding to S = S z = 0. Here S is the total 
spin of two electrons. The lowest energy excited states 
are the spin-triplet states with spatially antisym- 
metric wave- function corresponding to parity -1, and 

Here a = 0, ±1, and \T a ) is one of the spin-triplet states: 
|T +1 ) = I TT>, |T ) = (| U) + | lt»/\/2, and |T_ X ) = | ||) 
corresponding to 5 = 1, and S z = +1,0, — 1 respectively. 
Note that at B = (or zero Zccman coupling) these 
states are degenerate and thus ^ x is three-fold degener- 
ate. We define spin singlet-triplet energy gap of electrons 
A ( sl = E*^ - E { G ] where H e \^ G ) = £^ e) |* G ), and 

H e \*S>x) = Ex In tnis two level system there are 

two other excited states with spin-singlet ^f s ^(ri,f 2 ) — 
[Ptp+(ri)<p+(r 2 ) - aip-(fi)(p-(f 2 )] \S ), and 
*f (n,f 2 ) = 73 [^ + (fi) ¥ J_(f 2 ) + ^ + (f 2 ) ¥ J_(ri)] \Sq). 



trons perturbativcly. In this limit H em is assumed to be 
small, compared to unperturbed Hamiltonian H e . In the 
low magnetic field limit, the ground state of two electrons 
in coupled quantum dot is spin-singlet .It follows 
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Here G and X denote the ground and excited 
states of electrons in quantum dot systems, and 
X e {^x^x, 1 ^ ^ 1 }- The uniqueness (non- 
degeneracy) of the ground state has been assumed implic- 
itly. To obtain effective interaction between two Mn, we 
calculate the matrix elements of H em . It is straightfor- 
ward to show that (vf^ \H e 
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<&{Ri) = (p+iR^ip-iRi). We finally find 
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where A/// 



" 2A ( = 



Note that the Mn-Mn coupling for a lateral quantum dot 
molecule with two magnetic impurities localized at the 
center of each dot is given by 



A. Weak e-Mn coupling 



A12 = 



X\P sd tp 2 L {Ri)<p 2 R (R2) 
2AW 4(1 -W*) 



> 0. 



(7) 



Assuming the weak electron-Mn interaction limit, we 
calculate Mn-Mn effective Hamiltonian mediated by elec- 



Here we assume that -Ri and R 2 are the position of Mn's 
centered at left and right dots, and therefore the elec- 
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tron wave-functions at the opposite position of Mn's, 
viC^Ks) an d ^p 2 R {R\) are negligible due to high localiza- 
tion of the wave- functions. Because A12 is positive, the 
coupling between two Mn mediated by electrons is anti- 
ferromagnetic with M = as the ground state. For Mn, 
this state is separated by an energy gap, 30A 12 , from the 
ferromagnetic state M = Mi+M%(= 5). There are series 
of canted states with M = 1, ... ,4 between M = and 
M = 5. 



IV. EXACT DIAGONALIZATION 

The effect of V p at B = on both a and A^j, is shown 
in Fig. [T] This calculation is based on CI method us- 
ing 20 SP-orbitals (400 electronic configurations). With 
increasing V p the inter dot tunneling and the overlap be- 
tween L and R wave-functions decreases. This results 
to the decrease of a (hence A) and A^ simultaneously. 
The expansion of the ground state wave function |Wg), in 
terms of leading configurations of two electrons is shown. 
The contribution of the rest of configurations is negligi- 
ble. 

With increase of magnetic field, the electron spin 
singlet-triplet energy gap A^, decreases. Close to the 
transition point, Aoy vanishes, and the perturbation 
method fails. An unpertabative approach has to develop 
to calculate the low lying energy states of H in order to 
map H into H^ m . Here we exactly diagonalizc Hamil- 
tonian H by expanding the many body wave-function 
in the basis of electron-Mn configurations: |\&, M) = 
4rcL'l°) ® \M zU M z2 ). Because of [S z ,H em ] ± 0, 
(S is the total spin operator of two electrons), the 
states with different S z are mixed, hence the dimen- 
sion of matrix H that has to be diagonalized is given 
by N c = (2Mi + 1)(2M 2 + 1) £^ T=0 N S p[/[N^.(N SP - 
iV T )!]iVsp!/[^x!(iVsp - Ny)\\. N = N } + N± is the num- 
ber of electrons (here N = 2), and N$p is the number 
of single particle orbitals. To check the convergence of 
CI we perform exact diagonalization using single particle 
orbitals up to N$p — 20. The result of this calculation 
and the magnetic phase diagram of Mn is summarized in 
Fig. [2] where the electron spin singlet-triplet phase dia- 
gram is calculated by configuration interaction method 
for two electrons in lateral quantum dot molecules. 

Fig. [3] shows the lowest energy gap, A = Ej — E J=0 , 
calculated for two electrons and two Mn in lateral quan- 
tum dot molecule as a function of cyclotron frequency 
uj c = eB/m*c and J. Here J = M + S is the to- 
tal electron-Mn spin operator (M = M% + M%, and 
S = Si + S2 are total Mn and electron spin operators). 
For illustration we switched off the electron and Mn Zee- 
man couplings. Here the singlet-triplet transition occurs 
because of change in wave functions and e-e Coulomb 
matrix elements. Close to the transition point where the 
single particle energy levels of valence electrons are de- 
generate (half- filled) , the e-e Coulomb interaction leads 
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FIG. 3: (Color online) Energy gap of a system of two elec- 
trons and two Mn in lateral quantum dot molecule with 
V p = 7 as a function of cyclotron energy (magnetic field). At 
uj c = 1.55 electron spin singlet-triplet transition is seen. The 
vertical dashed line marks the transition point. This tran- 
sition induce para-ortho transition in Mn's. Below (above) 
this transition the ground state is identified by total angular 
momentum J = (J = 6). A magnified part of main figure 
is shown in inset. 



singlet-triplet transition in accordance with spin Hund's 
rule. The eigenvalues of H are grouped into J — 0, . . . , 6. 
States with given J are 2 J + 1-fold degenerate. It is con- 
venient to characterize these states based on the total 
spin of electrons, e.g., spin singlet (S — 0) and triplet 
(S = 1) and total spin of two Mn with M = 0, . . . , 5. In 
this work we are interested in the magnetic ordering of 
two Mn that can be described by anti-ferromagnetic, fer- 
romagnetic and canted states corresponding to M = 0, 
M = 5, and M = 1, . . . , 4. As it is shown in Fig. [3]spin of 
electrons undergo singlet-triplet transition at u>* = 1.55 
and V p = 7. Within cj c < w*, J = M = S = is 
the non-degenerate ground state. At lj c — lu*, the en- 
ergy gap of antiferromagnetic, ferromagnetic and canted 
states vanish all together and the ground state switches 
to ferromagnetic state with maximum spin multiplicity 
corresponding to J = 6, M = 5, and 5=1. In the limit 
of strong magnetic field there are 2 J + 1 = 13 degen- 
erate states that form the ground state. However this 
degeneracy is removed by Zeeman coupling that guaran- 
tees the uniqueness of the ground state with M z = —5 
and S z — —1. Within the resolution of our exact di- 
agonalization we did not observe any range of magnetic 
field that the ground state exhibits canted ordering. In 
Fig. [3] at B — and V p = 7 we compare the energy gap 
calculated using CI, Ej=q — Ej=o — 10 mRy*, with per- 
turbation approach. We find Agj, = 0.013, A12 = 0.38, 
hence Ej =6 — Ej =0 — 11.4 all in mRy* , in qualitative 
agreement with exact energy gap. 
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V. CONCLUSION 

In this work we studied phase diagram of quantum 
dot molecules consist of two electrons and two mag- 
netic impurities (Mn) confined in each dot. We demon- 
strated that the spin singlet-triplet transition of two 
electrons that are controlled by external electric gate 
voltage and magnetic field, can induce ferromagnctic- 
antifcrromagnetic transition in the magnetic impurity 
system. Therefore, Mn-Mn spin transitions mediated by 
e-Mn exchange interaction can be controlled indirectly 
by external electric gate voltage and magnetic field. This 
allows us to suggest application of spin of magnetic im- 
purities for the entanglement of qubits in quantum in- 
formation processing. The advantage of using spin of 
magnetic impurities as qubit, instead of QD electrons re- 
sides in the possibility in achieving higher spin coherence 
time. In analogous to the magnetic moment of nuclear 
impurities in host semiconductor systems, we speculated 
that the spin coherence time in the Mn system is ex- 
pected to be longer than the QD electron system due 
to small localization length of Mn d-orbitals that sup- 
presses the qubit spin-orbit coupling as well as hypcrfinc 
interaction with the magnetic moment of nuclei of host 
semiconductor. Our analysis based on exact diagonaliza- 



tion allows mapping the electron-electron and electron- 
Mn Hamiltonian H into an effective Mn-Mn Heisenbcrg 
Hamiltonian H^ m = Ai 2 Mi • M 2 in agreement with the 
perturbative results, e.g., an RKKY model calculated for 
weak coupling at low magnetic fields. Consistent with 
the magnetic field dependence of the lowest lying states 
of full Hamiltonian H , Ai 2 changes sign at critical mag- 
netic field that leads to spin singlet-triplet transition of 
two electrons in lateral quantum dot molecules. This is 
a level crossing that results to para-ortho transition in- 
duced by electrons in artificial H 2 molecules where the 
magnetic impurities resembling the magnetic moment of 
nucleus of the actual H 2 molecules, and the interaction 
between two Mn at high magnetic field is determined by 
the magnetic polaron effect. 
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